Partial decoherence in a qubit coupled to a control register subject to telegraph noise 
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A qubit (containing two quantum states, 1 and 2), is coupled to a control register (state 3), which 
is subject to telegraph noise. We study the time evolution of the density matrix p of an electron 
which starts in some coherent state on the qubit. At infinite time, p usually approaches the fully 
decoherent state, with p nm — 5 nm /3. However, when the Hamiltonian is symmetric under 1 2, the 
element p 12 approaches a non-zero real value, implying a partial coherence of the asymptotic state. 
The asymptotic density matrix depends only on Re[p 12 (t = 0)]. In several cases, the information 
stored on the qubit is protected from the noise. 
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Quantum computation operates on information stored 
in "qubits" , which are superpositions of two basic quan- 
tum states [l|. Clearly, quantum computation requires 
the stability of the quantum state stored on each qubit, 
and therefore it can be used only while this state remains 
coherent 0]. Interactions between qubits and their en- 
vironment, including input-output measurement devices, 
can cause decoherence which destroys the information 
stored in the qubits. Attempts to avoid decoherence have 
led to studies of decoherence-free subspaces, within which 
the quantum state remains protected 0]. The states in 
such subspaces are practically decoupled from the sources 
of the decoherence, due to symmetries of the system. 

Several proposals [1, Q to reduce the decoherence of 
the single qubit states use three basic states: two degen- 
erate states (|1) and |2)) forming the qubit, and the third 
state (| 3)) representing a control register. Models include 
an electron in a three-state atom [J], or in three single- 
level quantum dots @]. Below we refer to these states 
as sites in a tight binding model on a ring. The general 
Hamiltonian for such a system contains "site" energies 
e n and "hopping" matrix elements J nm , 

H Q = J2^n\n)(n\-(J nn+1 \n)(n + l\+h.c.)) , (1) 

n 

where we identify n = 4 with n = 1. In the non-stochastic 
case, one observes Rabi oscillations of the electron wave 
function, with frequencies equal to differences between 
the three eigen-energies of Tt . These oscillations are 
expected to decay when the system is coupled to the 
environment, due to decoherence. 

We represent the latter coupling by a stochastic tele- 
graph noise @, which affects only the control register: 
e 3 follows a Markov process or a two level jump process, 
and is given by e 3 = e — f(t)£, where f(t) jumps ran- 
domly from 1 to —1 (or from —1 to 1) with rate w . (or 

w, ). These jumps in f(t) result from a contact with 

some noise source, so that the probabilities p± to find 
/ = ±1 obey a Boltzmann distribution with energy gap 



AE at temperature T,p + = 1 —p_ = 1/[1 + e - AE / k B T ] 

Detailed balance implies p_w^ = p + w_ 

fore the jump rates can be written as 
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and thcre- 
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where A = +w , . Such a fluctuating site energy can 

be envisaged to occur because of noisy gate voltages or 
because of the coupling to a device (e.g. a point contact) 
which measures the occupation on that site Q. Being 
based on symmetry considerations, our results may be 
expected to hold in other models as well. 

Writing the quantum state of the system as \ip) = 
^2 n c n \n), the density matrix p = \ip}(4>\ has elements 
Pnm = c n c *m- Decoherence is usually associated with 
dephasing: the noise generates random fluctuations in 
the relative phases of c n and c m , causing the decay of 
the average off-diagonal matrix elements. Full decoher- 
ence corresponds to a fully "mixed" state, where all the 
off-diagonal elements of p vanish and all the three ba- 
sis states (in any representation) have equal occupation 
probabilities (in our case p nn = 1/3). In contrast, relax- 
ation describes the decay of the system towards a "pure" 
state, which is a coherent eigenstate (the ground state at 
zero temperature) of the non-stochastic Hamiltonian 0] . 
Possible measures of decoherence include 



V=l- Tr(p 2 ) 



(3) 



(one always has Trp = 1) or the von Neumann entropy, 



S = -Tr(plogp) 



(4) 



In these expressions, p represents a certain average over 
the stochastic noise (see below). Both measures vanish 
for a "pure" state, where the eigenvalues of p are equal 
to 1 and two zeroes, and have maxima (2/3 and log 3, 
respectively) for the fully decoherent state. Intermediate 
values of these measures require that off-diagonal ma- 
trix elements of p differ from zero, representing partial 
decoherence. Such deviations from full decoherence, as 
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found below, can probably be interpreted as a mixture of 
decoherence and relaxation 

Below we find the stochastic evolution in time of p, 
beginning with an initial coherent state in which the 
electron is on the qubit [i.e. R + = 1, where = 
Pn(t = 0) ± p 22 (t = 0)]- F° r arbitrary parameters of the 
non-stochastic Hamiltonian (e„, J nm ) and of the noise 
(£, !» T± ), our system develops full decoherence. How- 
ever, when the Hamiltonian is symmetric under the in- 
terchange of the two qubit states, 1 <-> 2, then the sys- 
tem maintains partial decoherence: although p 13 and p 23 
(i.e., the matrix elements which mix the qubit and the 
register states) approach zero, the matrix element p 12 
can approach a non-zero real value, implying an entan- 
glement of the final qubit states |1) and |2). Further- 
more, these states attain asymptotically equal popula- 
tions Pn = p 2 2 = (1 — P33)/2- Surprisingly, these asymp- 
totic values of p depend only on the initial value of the 
parameter y = Re[p 12 (i = 0)], which contains informa- 
tion on the initial entanglement of the two qubit states, 
and not on any of the other parameters (of the Hamil- 
tonian, the noise or the initial state)! A fully coherent 
initial state would have |p 12 | 2 = PnP 22 , anc ^ therefore 
y 2 < [1 — (i? _ ) 2 ]/4. In particular, the limiting values 
y = ±1/2 (and R~ = 0) correspond to the "bonding" and 
"anti-bonding" states (| 1) ± |2))/ -s/2, which are shown be- 
low to be protected from the noise, forming an example 
for a decoherence-free subspace. 

Here we discuss the simplest symmetric case, which 
contains only hopping between |1) or |2) to |3): we re- 
place H by H = H a - f(t)V, with V = C|3)(3|, and 
set e 1 = e 2 = 0, e 3 = e, J 12 = 0, J 13 = J 23 = J. 
The 3x3 density matrix p(t) obeys the von-Neumann- 
Liouvillc equation, id t p(t) — \H, p(t)]. It is convenient to 
rewrite this equation using the Liouvillc super-operators, 
id t p(t) = H x p{t) = [Ho - f{t)V x (t)]p{t), with the ob- 
vious definition A x p = [A, p(t)]. 

A treatment of the equation of motion with the under- 
lying stochasticity in f(t) can be found in the literature 
on the stochastic theory of lineshapes [lj| ■ Here we fol- 
low Blume [ll|, and replace each dynamic operator 0(t) 
by a 2 x 2 matrix 0(t), such that (6|0(i)|a) represents 
the conditional average of 0(i), given that the stochastic 
variable f(t) had the value a(= ±1) at time t = and 
the value b (=±1) at time t. At the end one may average 
over the stochastic process. The choice of the averaging 
procedure depends on the particular experiment. If one 
performs only one experiment, with a single initial state 
of the noise, then the initial value of f(t) is fixed (say at 
a). The final value of f(t) can then be either b = +1 or 
b = — 1, and thus an appropriate average is 

PUa = J2 ( b \ d \ a ) ■ (5) 
b=±l 

In many quantum measurements it is more appropriate 
to repeat the measurement many times, on systems which 



are prepared in the same way. In this case, it is necessary 
also to average over the initial value of f(t) (which equals 
a = ±1 with probability p a ), namely 

P}av = E Pa(b\0\a) EE J^PaPUa ■ (6) 
a,6=±l a 

For finite temperature T, i.e. \Ap\ < 1, where Ap = p + — 
p_, our calculations yield similar qualitative results for 
all these averages. In particular, all the averages coincide 
in the asymptotic limit t — > oo. 

We next construct a Markov chain of events in the time 
interval to t, where the back and forth jumps of f(t) are 
random and Poisson distributed [l2| . The time evolution 
of (b\p(t)\a) is now divided into two parts: 

d t (b\p\a) = - l [H x -bV x ](b\p\a) 

+ w b ,_ b {-b\p\a) - w_ b b (b\p\a) . (7) 

The first term on the RHS applies if f(t) remains un- 
changed at time t (i.e. equal to b). In this case, the time 
evolution proceeds with the Liouvillc operator which cor- 
responds to the original Hamiltonian, with e 3 = e — bC,. 
The last two terms arise if f(t) flips exactly at time t, ei- 
ther from — b to b (first term in the second line) or from b 
to —b (last term). Equation ([7]) can be put into a matrix 
form (in the stochastic 2— dimensional space), 

d t p(t) = H(W X T- V X F) + W]p(t) , (8) 

where 1 is the 2x2 unit matrix, while 

W = ( - W -+ W +~ )=X(f- T) (9) 

is the relaxation matrix [see Eq. In the above we 

have used the notations 

M°-°0- (10) 

Note that Eq. JS]) contains two types of operators: Hq 
and V x are Liouville super-operators (represented by 9 x 
9 matrices), while F and W are 2x2 matrices acting in 
the 2— dimensional space of the stochastic variable. The 
matrix p is of order 18 x 18. 

Using a Laplace transform, p(s) = / °° dte~ st p(t) (un- 
less specifically stated, all fPs below depend on s), Eq. 
([8]) becomes 

(si - W)(p n ± p 22 ) = i?±T + 2 JIm[p 13 ± p 23 ] , 

(si-W)p 12 =Pi 2 (0) -iJ(pi3 -P32) , 

[si - W - ie\p 13 = U (p 33 - p u -p l2 ) , 

[si- W - iS]p 23 = iJ(p 33 - p 22 - P21) , (11) 

where we have set p 13 (0) = p 23 (0) = and defined e 
f 1 — £F. Using the last two equations, and defining 

A= [sl-W-ie]- 1 = A R + iAj , (12) 
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we find Im[p 13 + p 23 ] = JA R B, where B = 2p 33 — p xx — 
p 22 — 2Rc[p 12 ]. Substitution into Eq. (jTTJ) then yields 

(si - W)(p n +p 22 )= R + l + 2J 2 A R B , 

(sl - W)Re[p 12 ] =yl + J 2 A R B . (13) 

It is easy to check that En Pnn CO = 1 implies that 

(Sl -W)J2 n Pnn = h ^- En Pnn = ^ whcrC ? = 

(sl — W)~ l is the Laplace transform of the probabilities 
to start at a and end up at b at time t. Combining this 
identity with the above equations, one has 



B=(2- 3R + - 2y) [sl - W + 8 J 2 A R }~ 



(14) 



To complete the solution of the Liouville equations we 
also hnd two coupled equations for p ±1 — p 22 and Im[/? 12 ]: 

Pu - P22 = C(R-\ + AJ 2 AjZlm[p 12 (0)}) , 

Im (Pl2] = ^( Im [Pl2(0)] - ^/(Pll - P22)) , 

Z = (si-W + 2J 2 A R )- 1 , 

C = (Z- 1 + U^AiZAiY 1 . (15) 

The time evolution of the density matrix follows in- 
verse Laplace transforms, which involve decaying oscil- 
lations, approaching asymptotic time independent limits 
as t — > 00. Before presenting an example of this time 
evolution, we consider these latter limits, using the iden- 
tity p(t — > 00) = lim ;! ^o[ s p(' s )]- Noting that Aj, Z 
and C all remain finite as s — ► 0, we conclude that 
p u (oo) - p 22 (oo) = 3[p 12 (oo)] = 0, p 13 (oo) = p 23 (oo) = 
0. Tedious but straightforward algebra is needed to show 
that 



lim J 2 A R [sl -W 

S— >0 



8J 2 A R Y 



1 

Tg 



l + e/C 1 + e/C 
1-e/C 1-e/C 



(16) 



Using also P(t — > 00) = lim s _»o[s-P] = T, we find that 
the matrix p splits as 



( b \p mn \a) = p™ n (b\T\a) = p™ n p b 



(17) 



with 



R+ + 2-2y 6y + 2-3R+ 
6y + 2-3R+ R+ + 2-2y 

4 - 2R+ + Ay 

(18) 



Since the factor p b in (JTTJ) gives the probability to find 
the stochastic variable in state b, it does not affect the 
actually measured density matrix, which will always ap- 
proach p°°. In particular, [T] av a = 1 and therefore all 
the averages will yield the same results. When the elec- 
tron is initially placed on the qubit, i.e. R + = 1, then 



p°° depends only on y, and not on the other details of 
the initial state. The result is also independent of the 
non-stochastic Hamiltonian (e, J) and of the stochastic 
noise (Ap, £, A), as long as A > and |Ap| < 1. In 
contrast, when Ap = ±1, i.e. at zero temperature, the 
system jumps to the state with e 3 = e =p £ and then it 
stays there forever, behaving like the non-stochastic sys- 
tem with that energy. 

Substituting Eq. (fl~8|) with R + = 1 into the decoher- 
ence measures ([3]) and (J4]) yields maximal decohcrence, 
p nm (oo) = S nm /3, when y = 1/6. When the electron 
starts in state |1), then y = 0, and one has asymp- 
totic partial decoherence, p n = p 22 = 3/8, p 12 = —1/8. 
When the initial state of the qubit is fully entangled, 
(|1) +e ra |2))/V2, then 



P(0) = \ 



1 







(19) 



so that y = icosa and = i cos 2 (a/2). Measuring 
this occupation of the control register then yields a mea- 
sure of the phase a, in spite of the decoherence! 

Both measures of decoherence, Eqs. (O and ([1]), vanish 
at all times if the electron starts with y = —1/2 (a = tt), 
which corresponds to the "anti-bonding" state \AB) = 
(|1) — |2))/v2- Being fully antisymmetric under 1 *-> 
2, \AB) does not couple to |3), and is not affected by 
the symmetric Hamiltonian. Therefore, this state forms 
the decoherent free subspace, which remains constant and 
fully protected at all times. As y increases from —1/2, 
one encounters a mixture of this "anti-bonding" state and 
the other two orthogonal eigenstates of Ho, and thus the 
system develops various degrees of partial decoherence. 
The value y = 1/2 (a = 0) corresponds to the initial 
"bonding" state, \B) = (|1) + |2))/V2. In this case, pf{ = 
P22 = Pu = 1/4) and the projection of the quantum state 
onto the qubit subspace remains in the fully coherent 
state \B). 

To demonstrate typical results, we start with an arbi- 
tr ary choi ce of parameters, y = —0.3, R + = 1, R~ = 
\/l — Ay 2 and set e = J. Figure [TJ shows the three occu- 
pations p nn (t). The non-stochastic limit (left panel) ex- 
hibits Rabi oscillations with period 2ir/J, corresponding 
to the difference between eigen-energies of TC ; the cou- 
pling of the qubit to the control register causes periodic 
partial periodic occupation of the latter. In the stochas- 
tic cases, small values of A(<C J) and/or £(<C J) yield 
small corrections to this non-stochastic picture at short 
times, with an eventual "flow" to the asymptotic behav- 
ior described above. The remaining panels of Fig. [JJshow 
[p] av for a relatively large stochasticity, £ = A = J/2 and 
Ap = p + -p_=. 0.1 (high T) or 0.9 (low T). The other 
averages [p] av ± exhibit similar qualitative behavior, and 
all the averages approach the asymptotic values given by 
Eq. (fT5|) . Clearly, the decoherence time is longer for the 
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lower temperature (larger |Ap|). Figure [5] shows the de- 
coherence measure T>, calculated using the three averages 
\p\av ± (dashed lines) and [p] av (full lines). Note also that 
the asymptotic limit T> = 0.34 is much smaller than the 



fully decoherent value 2/3, demonstrating the partial de- 
coherence. Also, the approach to this asymptotic value 
is much slower for the lower temperature. 
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FIG. 1: The site occupations for an initial coherent state with y — —0.3, and e = J. Left: non-stochastic limit. Center and 
right: average occupations for stochastic cases, with f = A = 0.5 J, Ap = .1 and .9. pn and p 2 2 approach (3 — 2j/)/8 = 0.45 
and p 3 3 approaches (1 + 2y)/4 = 0.1 (shown by horizontal lines). 
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FIG. 2: Time dependence of the decoherence measure T>, for 
the same parameters as in Fig. [1] Left (right): Ap = 0.1 (0.9) 
and with f = A = 0.5J. Small (large) dashes correspond to 
[p]au,+ ([p]av,-)- Full lines correspond to [p] av . The horizontal 
line shows the asymptotic value, T> = 0.34. 

The equations of motion become more complicated 
when one adds direct hopping between the two qubit 
states, i.e. J 12 ^ 0. We have solved the equations of mo- 
tion for this generalized case, by inverting the full 9x9 
matrices for the super-operators [l3j]. Surprisingly, the 
asymptotic result (|18|) remains the same, independent of 
J 12 . This robust result depends only on the symmetry 
1 <-> 2. We have also checked that any deviation from this 
symmetry (e.g. by introducing a Aharonov-Bohm flux 
inside the ring formed by the three sites, or by having 
different energies on sites 1 and 2) immediately changes 
the asymptotic limit back into the fully decoherent one. 

It would be very interesting to check if Eq. (fl8|) also 
holds for other sources of noise, e.g. a capacitive coupling 
to a point contact Q. It would also be interesting to 
consider generalizations of these results to systems con- 
taining several qubits. We hope that our results will also 
be tested experimentally 
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